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Abstract 

The secondary spectra of the gravitons induced by a waterfall-like field are computed and 
the general bounds on the spectral energy density of the tensor modes of the geometry are 
translated into explicit constraints on the amplitude and slope of the waterfall spectrum. The 
obtained results are compared with the primary gravitational wave spectra of the concor- 
dance model and of its neighboring extensions as well as with the direct Ligo/ Virgo bounds 
on stochastic backgrounds of relic gravitons. Space-borne interferometers (such as Lisa, Bbo, 
Decigo) seem to be less relevant but their potential implications are briefly outlined. 
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1 Motivations 



There are situations, in the early Universe, where waterfall-like fields arising in hybrid infla- 
tion [T] are amplified with spectral slopes which are even steeper than the ones characterizing 
the vacuum fluctuations. In the latter case the inhomogeneities of the scalar field induce a 
secondary relic graviton spectrum as opposed to the primary spectra induced directly though 
the coupling of the tensor modes of the geometry to the the space-time curvature. Provided 
the slopes of the secondary spectra are sufficiently steep, they can dominate the energy 
density of the produced gravitons at high frequencies. In this paper the secondary tensor 
spectra will be computed within a sufficiently general parametrization so that the obtained 
results can be of wider applicability. Recalling that the two point function of a canonically 
normalized scalar field is customarily defined as 

(a(k,r)a(p,r)) = ^-V a (k,r)5^(k + p), (1.1) 

the power spectrum of a can be parametrized as 

— ) T(kr), (1.2) 

"'max / 

where A a has the dimension of an inverse length and n a is the spectral slope. In Eq. (11. 2p the 
dimensionless function ^(kr) accounts for the time dependence and the power specctrum is 
normalized at the comoving wavenumber /c max corresponding to a comoving frequency ^ max 

Z/max = ^r = - 417/3 WJ U.28 x 10-0 U.15X10-0 (L3) 

where the fiducial valued are the ones corresponding to the WMAP 7yr data alone [2] 
supplemented by a tensor component with rx = 16 e < 0.36 (with e = —H/H 2 the standard 
slow- roll parameter); the parameter (3 appearing in Eq. (11.31) depends upon the width of the 
transition between the inflationary phase and the subsequent radiation dominated phase. In 
the notations of Eqs. (II. ip and (II. 2p the scale-invariant limit corresponds to n a = 1. The case 
n a = 3 characterizes the slope of quantum (vacuum) fluctuations. If n a > 3 the spectral slope 
is steeper than in the case of vacuum fluctuations. The amplified spectrum characterizing 
the waterfall field in hybrid inflation leads, according to recent analyses [3], to n a c± 4. The 
curvature perturbations at the comoving scale k p = 0.002 Mpc -1 are characterized by an 
amplitude An 

(n(k,r)n(p,r)) = —V 1l (k)5^(k+p), v n {k)=An[ Y ) • (1.4) 

In the notation of Eq. (I1.4p the curvature perturbations induced by the waterfall field can 
be estimated for scales comparable with k p and the results are such that An — O(10 -54 ) 

2 In Eq. (|1.3j) Ar denotes the amplitude of the curvature perturbations (see also Eq. (|1.4jl ): VLqx denotes 
the present value of the energy (or matter) density parameter in critical units and for the species X. 
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[3]. The latter figure could be compared, for instance, with the amplitude of the curvature 
perturbations induced by the adiabatic mode are An = (2.43 ± 0.11) x 10~ 10 with n s = 
0.963 ± 0.014 according to the WMAP 7yr data alone analyzed in the light of the ACDM 
paradigm and without tensors [2]. 

The smallness of the induced curvature perturbations at the scale k p (corresponding to 
a typical multipole i v ~ 30) guarantees the absence of modifications in the temperature 
and polarization anisotropies for all the multipoles relevant in CMB studies. The smallness 
of the scalar fluctuations of the geometry in the frequency range probed by CMB physics 
[3] does not exclude the possibility of a sharp increase of the spectral energy density of 
the tensor modes of the geometry at higher frequencies. Indeed, as it is well known, the 
typical frequencies probed by CMB physics correspond to low-frequency gravitons in the 
aHz range while i/ max is rather in the GHz rang^f). In the concordance model the tensor 
modes are, strictly speaking, absent. It is however possible to include them by adding, in the 
simplest situation, only one parameter (i.e. r T , already introduced after Eq. (11. 3p ) which 
controls, simultaneously, the tensor spectral amplitude and the tensor spectral slope. For 
frequencies in the GHz range and in the minimal tensor extension of the ACDM model the 
spectral energy density of the relic gravitons can be easily 0(1O~ 16 ) in units of the critical 
energy density and for tt — O(10~ 3 ). Can this figure change (even dramatically) when the 
secondary contribution to the relic graviton background is consistently taken into account? 
How large can the combined spectral energy density be at high frequencies? How does the 
resulting graviton spectral energy density compare with the other potential signals arising in 
slightly different extensions of the ACDM paradigm? These are few of the questions which 
will be specifically addressed hereunder. The purpose of section [2] is therefore to compute 
the slope and amplitude of the secondary graviton spectra. In section [3] the results of [2] will 
be used to ascertain wether and in which region of the parameter space the current bounds 
on the stochastic backgrounds of relic gravitons are adequately satisfied. In section 0] the 
secondary graviton spectra are compared with the primary signal induced directly by the 
amplified quantum fluctuations of the tensor modes of the geometry. Section [5] contains the 
concluding considerations. 



2 Secondary spectra of relic gravitons 

In hybrid inflationary models pQ the scalar potential can be written in terms of the slowly 
rolling inflaton field <p and in terms of the waterfall field a endowed with a mexican hat 
potential and directly coupled to the inflaton: 

W(<p, a) = - A (2ml - v^a 2 ) 2 + -f<f + ^V. (2.1) 
3 Recall that 1 aHz = l(T 18 Hz while 1 GHz = 10 9 Hz. 

4 The concordance model will be taken as a synonym of the ACDM scenario where A stands for the dark 
energy component and CDM for the cold dark matter contribution. 
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The true vacuum of a arises for a 2 = 2m^./vA and <p — while the false vacuum occurs for 
a = 0. The spectrum of a close to the background solution cr = falls into the class described 
by Eqs. (II. ip and (II. 2p [3J. At <p\ = 2y/\m 2 /g 2 the waterfall field becomes almost suddenly 
massless, and the (p)reheating phase start. To keep the discussion reasonably general, the 
case of a waterfall-like degree of freedom characterized by a steep power spectrurrjf] n a > 2 
will be considered. The latter assumptions are consistent with the sudden approximation 
where the transition between the inflationary and the radiation-dominated phase is almost 
instantaneous. In the case of the sudden approximation the parameter of Eq. (II. 3p is 
exactly equal to 1 but it can become different from 1 if the transition to radiation is not 
sudden as computed in explicit numerical examples [I] . The tachyonic instability leading to 
the graviton spectra discussed in the present paper might also arise, under certain conditions, 
in inflationary models motivated by supergravity In this sense waterfall-like field (not 
necessarily arising in standard hybrid models) could be touched by some of the present 
considerations which will be kept, in what follows, as general as possible. To compute the 
spectral energy density of the relic gravitons the relevant quantity is the anisotropic stress 
of the field a, i.e. 

n(f , r) = ditrdja - ^5 i3 (d k a) 2 , (2.2) 

whose projection on the two tensor polarizations determines the spectral energy density 
of the relic gravitons. Consider, therefore, a conformally flat Friedmann-Robertson- Walker 
(FRW) geometry supplemented by its tensor fluctuations 5 t gij, i.e. 

V = a V)^ ^9ij(r,q) = -a 2 (T)hij(q,T). (2.3) 
The two polarizations of hij(q,r) obey, respectively, 

hi + 2Ub! @ + = 24n®(g, r), (2.4) 
h% + 2Hhl + q 2 h® = 2£ 2 Mq, r), (2.5) 

where the prime denotes a derivation with respect to the conformal time coordinate r and 
"H = a' I a. In Eqs. (12. 4 p and (12. 5p the Planck length is defined as 

£ p = = 1/Mp, (2.6) 

while, always, in Eqs. (12 .4p and (12. 5p II e (g, r) and n (g, r) denote the projections of the 
anisotropic stress on the standard basis of the tensor polarizations, i.e. 

n®(g, r) = 2^3/2 / ^ k e e(<?) k i k j r ) ~ £ r )> ( 2 - 7 ) 

n®(g, t) = 2(2^)3/2 / d * k e ®($ ki k i a ( k > r ) a (v~ k > r )' ( 2 ' 8 ) 

5 Even if the case n a > 2 is the one directly related to the waterfall spectra, the value of the spectral 
index will be kept generic to allow for a wider applicability of the obtained results. 
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where e®(g) and efjq) are defined as§ 

e ?j(<l) = - bibj), efjiq) = (a*&j + kdj). (2.9) 

As already mentioned, it will be assumed that the relaxation of a takes place right at the 
end of inflation; later on the background geometry is dominated first by radiation and then 
by matter. Over the frequencies affected by the radiation-matter transition the leading 
contribution to the spectral energy density comes from the amplified vacuum fluctuations 
which left the Hubble radius during inflation and re-entered during the matter-dominated 
epoch (see section H]). Conversely the leading contribution to the hight-frequency part of the 
spectrum can be computed from Eqs. (12. 4 p and f 12 . 5 j) . More precisely, during the radiation- 
dominated stage of expansion each polarization obeys the equation 

2£ 2 r 

h(xMr) = J o a(£)Il w (q,£)G(q,T,£)d£, (2.10) 

where A = ©, ® stands for each of the two polarizations and G(q,r,£) = sin [q{r — £)] in 
the case of a radiation-dominated Universe. To compute the spectral energy density of the 
relic gravitons we need an explicit expression for the energy density. In a FRW background 
of Eq. A2.3H the action of the tensor modes can be written as (see, for instance, [4]) 

Sgw = g^f / (PxdTy/^jjg^d^dpti*. (2.11) 

The energy-momentum tensor derived from Eq. (12. lip by varying with respect to g a/3 , as 
argued by Ford and Parker [6], is the correct quantity to use for the calculation of the spectral 
energy density: 

l l 

(2.12) 



V = — 



Alternatively one could use the energy-momentum pseudo-tensor j7] appropriately general- 
ized to curved backgrounds [H |9]. The energy-momentum pseudo-tensor defined from the 
second-order variation of the Einstein tensor and the expression of Eq. (I2.12p are formally 
slightly different and their peruse in the context of stochastic backgrounds of relic gravitons 
has been carefully gauged in JTUJ [11] (see also jl]): the two computational strategies lead 
exactly to the same spectral energy density for wavelengths shorter than the Hubble radius 
at each corresponding epoch (i.e. qr 1) and to quantitatively compatible results (within a 
factor 2) in the opposite limit. For wavelengths shorter than the Hubble radius the relation 
between the spectral energy density (in critical units) and the power spectrum is given by0 

^Gw(s,r) = = -^-P T (g,r), (2.13) 

p crit ding 12H 2 



6 Note that a, b and q form a triplet of mutually orthogonal unit vectors. 

7 In the present paper In denotes the natural logarithm while log denotes the common logarithm. 
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where p crit = 3H 2 /(a 2 £p) is the critical energy density and where P^i^q.r) is the power 
spectrum satisfying 

(Mx,t)/^(£ + f,r)> = J d\nqP T (q,r) S -^-. (2.14) 

The direction of propagation of the graviton q can be fixed along thez axis; consequently 
denoting with /i = cos?? the explicit form of Qcw(q, r) can be written, from Eq. (12.101) . as 

"swfex) = ^(^) 4 t~**j\ 4.(1 - ^.W^^A^r) (2.15) 



where 



p/ -'fcp "/-l |g — fc| f 



J (1) (fc, g, r) = f f ^ sin [g(r - 0] sin (*£) sin [|g - fc|£]. (2.16) 
Equation (12. 16|) follows by inserting, into Eqs. f !2.Tj) and (12. 8p . 

— * 

<r(A;, r)a(q- k, r) = a(k)a(q- k) — ^ — tlz___U.. (2.17) 

A;r |?-Jfe|T 

Denoting with p the combination \q — k\ the integral of Eq. fl 2 . 1 6 [) can be expressed as 

l (1) (k,q,r) = S ^MS [ Ci[(k- P + q)T}-Ci[( q - P -k)T} 

+ Ci[(q + p - k)r] - Ci[{k + V + q)r] + In [ ^Z^y ] } 

+ ^M{&[{k+p + q)T]-Si[(p + q-k)T] 

- 8i[(k -p + q)r] -Sip+p-g)r]|, (2.18) 



where, with standard notations [T2| [13]. 

ci( z) = - r ^i d t, si(z) = r —dt. (2.i9) 

J z t J z t 

The small and large argument limit of Ci(z) and Si(z) is given, respectively, by 
Ci(z) = 7 + In z + 0{z 2 ), Si(z) = ^ + C>(z 3 ), 

_ . . SU1Z COS Z , 7T COS 2 

Ci(z = + O ^ 3 , Si * =tt + ^ (z ~ 2 . 2.20 

z z z 2 z 

Note that z ^> 1 all the wavelengths are shorter than the Hubble radius and, in this limit 

7T singr f Akq 2 p 

Mi){k,q,T) = --cos qr + 



+gln 



4 * 4(gr)l [(q + k) 2 - p 2 ] [(q - k) 2 - p 2 } 
r q 2 - (k + p) 



q 2 — (k — P y 



The leading contribution to Eq. (I2.15P is given by the first term (proportional to cosgr) 
while the other terms are suppressed as qr > 1. The oscillations inside the Hubble radius are 
typically averaged in the final result by replacing cos 2 qr — > 1/2 (see, e.g. [HI 03]). Following 
the same procedure, the final expression for the spectral energy density of the relic gravitons 
is given by: 

■2 



^gw(<7,t) 



7T 



A n 



3 r k 



k 



Z(2){q/qm a *,k/k max ,n a ) d(k/k B 



192 \MpJ \q mi 
where, defining x = k/k m3uX and y = q/q max , 

X {2) (x,y, n<7 ) = j\l-i2 2 ) 2 [y 2 + x 2 -2x W } ( - n °- (i)/2 dii. 



(2.22) 



(2.23) 



The integrals of Eq. (I2.22p and (12.231) can be performed numerically and, in some approxi- 
mation, also analytically; the final result can be written as 

r 2 



^gw(<7,t) 



7T 



/ Ag 



180 VMt 



Q 



(In 



2(n CT -l) 



6 - n n 



+ 



(n a + 1)(5 - 2n CT ) 2n a -5\q- 



Q 



5— 2n„ 



. (2.24) 



Let us now pause for a moment and discuss the approximations used to derive Eq. (I2.24j) . 
The first remark is that the integral of Eq. (I2.23P seems to lead to a logarithmic divergence 
in the case n a = 4. Even in the case n a = 4, Eq. (I2.24p is in very good quantitative 
agreement with the numerical result. The rationale for this occurrence is that Eq. ( 12.221) 
can be explicitly written, in the case n a = 4, as 



G{y) 



192 \M P , 

1 f 1 ^{ 3(x 2 _ y2) 4 _ Axy{x 2 + y2)(3x 4 _ Ux 2 y 2 + 3y 4 } 

x + y\ 



96y 



+ 6(x 2 -y 2 ) A \n 



\x - y\ 



(2.25) 



(2.26) 



and where y p = k p /k m3uX = z/ p /z/ max .The integral appearing in Q{y) can be easily performed 
numerically by recalling that 



y p = 7.41 x 10 



-27 



1 



/3V0.01 



-1/4 



2.28 x 10" 9 



0"R0 



4.15 x 10" 5 



(2.27) 



To obtain a quantitatively accurate expression for the integral it is useful to write G(y) as 



G(y) 



96 



F(u) + [ 



1 dxF(s) 



y X S c 



rv dx 
'y P x 

where u = x/y and s = y/x while for a generic variable t the function F(t) is defined as 

II - t 



2.28) 



F(t) = 4t(l + t 2 )[Ut 2 - 3(1 + t 4 )] - Q(t 2 - l) 4 ln 

7 



1 + t 



(2.29) 



Since u < 1 for y p < x < y and s<lfory<a:<l, G(y) can be evaluated by expanding F(t) 
in the two domains and by keeping the first few terms of the expansion and by integrating 
the obtained result. In the limit y p — > the final result for Qgw(9; t ) coincides with the 
expression given in Eq. (12.241) for n a = 4. 

The approximation scheme leading to the result of Eq. (I2.24p is rather accurate both 
for the present purposes and in more general terms. To scrutinize more carefully the latter 
statement the approximate results based on Eq. (I2.24p will be compared with the explicit 
numerical integration in various situations. For this purpose it is appropriate to generalize 
G(y) to the case where n a ^ 4: 

G(y,n a )=y 3 C x n °l (2) (x, y, n a ), (2.30) 
J y P 

which reduces to Eq. (|2.26[) when n a = 4 as it can be easily verified with some algebra. The 
integral of Eq. (12.301) is simply regularized in the ultraviolet by selecting a cut-off and this is 
the rationale for the upper limit of integration. The results of the comparison between the 




log y log y 



Figure 1: The comparison between the numeric and the analytic results for the calculation 
of G(y,n a ). 

numeric evaluation of G(y, n a ) and the analytic expression of Eq. ( I2.24j) is illustrated in Fig. 
CD In the plot at the left the stars, the squares and the naughts reproduce, respectively, the 
numeric results for the different values of n a . The lines reproduce instead the analytic result 
obtained on the basis of Eq. ( 12.24]) . In the plot at the right the results for the case n a = 4 
are shown in greater detail for y close to y — 1. The full line of the right plot shows the 
numeric result while the dashed and dot-dashed lines illustrate instead the analytic estimate. 
More specifically the dashed line shows exactly the result of Eq. (12.241) while the dot-dashed 
line is based on the improved approximation where the next-to-leading term is kept when 
expanding F(u) and F(s) in Eq. (12.281) . i.e. 



G(y,n a ) = f/ {n °- 1] { 



n a -Q y 



5—2n a 



(2n CT -5)(n CT + 1) 2n a 



[n a - l){n a - 6) 
14 

The first line in Eq. (12. 3 1 p leads directly to the result of Eq. ( I2.24p while the second 
line represent the correction. The goodness of the agreement between the analytic and the 
numeric results is therefore quantitatively evident from Fig. [TJ Every approximation can 
be improved around the transition frequency k max by considering the possible exponential 
suppression of the modes k > k max as implied, in the case of gravitational particle production, 
by the so-called adiabatic theorem [16]. Similar agreement is reached when the cut-off is 
replaced by the exponential tail. It is however important to notice (as noticed in a related 
context [3]) that the width of the transition between inflation and radiation affects the 
definition of z/ max itself so that 1//3 increases from 1 to 6.33 in the example of [3]. These 
kinds of inaccuracies imply that the exponential tail is itself necessarily inaccurate and must 
be determined in the context of a specific model. Put in other words it is always possible 
to produce more accurate estimates by purportedly invoking a more accurate description of 
the transition regime. Since in this (as well as in similar cases) the understanding of the 
physics is approximate, we shall be content of the present approximation scheme which has 
the advantage of being not only accurate but also transparent. 



n a - 1U 



+ 



(n„ + 3)(2n„-7) 2n„ - 7 



(2.31) 



3 Observational constraints 

The results of Eq. (I2.24p estimate the secondary contribution to the spectral energy density 
of the relic gravitons whose primary signal emerge because of the direct coupling of the 
tensor modes of the geometry to the space-time curvature. The latter coupling amplifies an 
initial vacuum fluctuation during the inflationary quasi-de Sitter stage of expansion and the 
resulting spectral energy density can be written as 



hlQ GW (u,T ) = A/" T T^(z//z/ eq )r T ( — ) e 2 -^, (3.1) 

' h 2 Q R0 \ / A n \ 
.4.15 x lO" 5 / V2.28 x lO" 9 /' 



(It 



At = 3J4xlO- 15 f-^% =? V— ^^W), (3.2) 



where rx and ut are defined as 

r T = ^-y, n T = - 2e + T ln(,/, p ), a T = - 



[n s - 1) + ? 



(3.3) 



If qt = the tensor spectral index does not depend upon the frequency and this is 
the case which is, somehow, endorsed when introducing gravitational waves in the minimal 
tensor extension of the ACDM paradigm [2]. If ^ the spectral index is said to be 
running at a rate that depends upon the value of the scalar spectral index n s . Note that, by 
definition, rx represents the ratio between the tensor and the scalar power spectrum at the 
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scale k p . In Eq. ( 13.1 p Tn(v/v eq ) denotes the transfer function of the spectral energy density: 



Tn{vM = )jl + c L \^)+b 1 (j^) , c 1 = 0.5238, b x = 0.3537, (3.4) 
where 

= = 1.288 x W^mo) ( h ° Q ™ Y 1/2 Hz. (3.5) 

q 2vr V0.1334/ V4.15 x 10- 5 / V ; 

Equation H3.4[) has been derived in [10] (see also [U [11]) but other authors prefer to work 

with the transfer function of the power spectrum (see, e. g. [14j HE]). The two procedures are 

equivalent but the transfer function of the spectral energy density is easier to assess directly 

in numerical terms by integrating the evolution of the mode functions across the matter- 

ratiation transition for a given initial spectrum. Conversely if one oughts to evaluate the 

transfer function for the power spectrum Pt(v, r), the obtained result must then be connected 

to Oqw 

(z/, r) and this should be done carefully since the power spectrum oscillates and the 
oscillations must then be appropriately averaged [I] . The specific values of the cosmological 
parameters determined using the WMAP 7yr data alone in the light of the vanilla ACDM 
scenario are: 

(fi b , fi c ,fi de , h , n s , e re ) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088). (3.6) 

If a tensor component is allowed in the analysis of the WMAP 7yr data alone the relevant 
cosmological parameters are determined to be 

(fi b , fi c ,fi de , ho, n s , e re ) = (0.0430, 0.200, 0.757, 0.735, 0.982, 0.091). (3.7) 

In the case of Eq. ( 13. 6p the amplitude of the scalar modes is An = (2.43 ± 0.11) x 10~ 9 
while in the case of Eq. ( 13. 7p the corresponding values of An and of r T are given by 

An = (2.28 ±0.15) x 10~ 9 , r T < 0.36 (3.8) 

to 95 % confidence level. To avoid confusions it is appropriate to spend a word of care on 
the figures implied by Eqs. (13.71) and (13.81) which have been used in the numeric analysis just 
for sake of accuracy. At the same time the qualitative features of the effects discussed here 
do not change if, for instance, one would endorse the parameters drawn from the minimal 
tensor extension of the ACDM paradigm and compared not to the WMAP 7yr data release 
but rather with the WMAP 3yr data release, implying, for instance, An = 2.1^2 3 x 10~ 9 , 
n s = 0.984 and rx < 0.65 (95 % confidence level). With the previous caveats the spectral 
energy density can be written, in explicit terms, as 



hlQ GW (u,T ) = T^(u/u ccl )lAf T r T ( — 

+ K 



2- 



g "max 



p ' 

2(n CT -l) r g _ n ^ I ( V \ 5 ~ 2n ° 



n a + 1)(5 - 2n a ) 2n a - 5 



(3.9) 
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where 

It is interesting to remark that, for a given (fixed) ratio u/u max 



lim 

, CT ->5/2 



?v — 6 1 / v 



(n CT + l)(2n CT -5) 2^-5 



5— 2n CT -i 



2 

- + In (z/ max /V), (3.11) 



so that the apparent divergence occurring for n a = 5/2 in Eq. (13. 9 p just lead to a logarithmic 
enhancement. The same kind of limit arises in the next-to- leading contribution appearing 
in the second line of Eq. ( I2.3ip . i.e. 



lim 



n n — 10 1 ( v 



{n a + 3) (2n a - 7) 2n CT - 7 \u, 



7-2n a 



— + ln(iw/i/). (3.12) 



Denoting with £ the order of the perturbative correction to the leading result of Eq. (13. lip , 
the logarithmic enhancement arise for n a = (5 + 2£)/2 and their magnitude is 



lim 

>(5+2C)/2 



n a - 6 - 4C 



(n CT + l + 2C)(2rv-5-2C) 



2n a -5-2CVz/ 1 



7/ 



5+2C-2n 



2 



7 + 6C 



+ ln(z/ max /z/). (3.13) 



The relic graviton spectra must satisfy various constraints which will also imply, in par- 
ticular, quantitative bounds on the parameters of the secondary spectra. The spectral energy 
density of the relic gravitons must be compatible not only with the CMB constraints (bound- 
ing, from above, the value of r<r) but also with the pulsar timing bound [18j EE] and the 
big-bang nucleosynthesis constraints [271 EHl EH] . The result of Eq. ( 13.1 ft automatically in- 
corporates the constraint stemming from the CMB observations encoded in the value of pp. 
As we move towards larger frequencies the pulsar timing bound demands 

fi(^puisar,r ) < 1.9 x 10" 8 , z/ pulsar ~ 10 nHz, (3.14) 

where v VVl \ S w roughly corresponds to the inverse of the observation time during which the 
pulsars timing has been monitored. 

Following the analyses on the stochastic relic graviton backgrounds [20] (see also [21]) let 
us define z/lv = 100Hz as the Ligo/Virgo frequency. In Ref. [22] (see also J23J, [23] ) the current 
limits on the presence of an isotropic background of relic gravitons have been illustrated. 
According to the Ligo collaboration (see Eq. (19) of Ref. [23]) the spectral energy density 
of a putative (isotropic) background of relic gravitons can be parametrized as: 



i a+3 

V 100 Hz, 



Oqw (y, T ) = GW a ( " ) (3.15) 



For the scale-invariant case (i.e. a = — 3 in eq. ( I3.15P ) the Ligo collaboration sets a 90% 
upper limit of 1.20 x 10~ 4 on the amplitude appearing in Eq. ( I3.15p . i.e. ficw-3' Using 
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different sets of data (see [221 121]) the Ligo collaboration manages to improve the bound 
even by a factor 2 getting down to 6.5 x 1CT 5 . Recently this result has been improved to 
6.9 x 10~ 6 [26] and this is the figure used in the forthcoming plots when generically referring 
to the Ligo/ Virgo bound. 

Assuming that there are some additional relativistic degrees of freedom, which also have 
decoupled by the time of electron-positron annihilation, or just some additional component 
p x to the energy density with a radiation-like equation of state (i.e. px = Px/3), the 
effect on the expansion rate will be the same as that of having some (perhaps a fractional 
number of) additional neutrino species. Before electron-positron annihilation we have px = 
(7/8)AN u p 1 and after electron-positron annihilation we have px = (7/8)(4/ll) 4 / 3 AN U p 7 ~ 
0.227 AN U p y . The critical fraction of CMB photons can be directly computed from the value 
of the CMB temperature and it is given by /iq0 7 = p 7 / p crit = 2.47 x 1CT 5 . If the extra energy 
density component has stayed radiation-like until today, its ratio to the critical density, Qx, 
is given by 

h&x = hfj- = 5.61 x 10- 6 A^( 2 ^ _ 5 ). (3.16) 
If the additional species are relic gravitons, then [27] |28j 129] : 

h 2 £ m Jtt GW (v,T )d\nv = 5.61 x 10 ~ 6A ^( 24 7°^° -5 ). ( 3 - 17 ) 
where i/ max is given by Eq. (11. 3p and z/bbn is given by: 

u hbn = 2.252 x 10- 11 ( V4 ( ( k2 ° QR0 - ) V4 Hz ~ 0.01 nHz, (3.18) 
U0.75/ V MeV/ V4.15 x 10- 5 / ' V ; 

In Eq. (I3.18P g p denotes the effective number of relativistic degrees of freedom entering the 
total energy density of the plasma. While z/ eq is still close to the aHz, i/bbn is rather in the 
nHz range. In Fig. [2] the limits on the amplitude A a have been illustrated as they arise 
from Eq. (I3.17P and for different values of AN V . In the plot at the left of Fig. [2] the two 
parameters are A a and n a while in the plot at the right, on the horizontal axis, parametrizes 
the indetermination on the maximal frequency or, in a complementary perspective, the 
limitations of the sudden approximation. In Fig. [3] the current Ligo/Virgo bounds are 
contrasted with the expected sensitivities of some of the planned space-borne interferometers 
such as Lisa [30] and the Bbo/Decigo [3T] projects. In the Lisa case the plots have been 
obtained by assuming a nominal sensitivity of 10 -11 (for the spectral energy density) and 
in a frequency range centered around the mHz, i.e. i>usa = 10~ 3 Hz. In the Bbo/Decigo 
case the sensitivity has been taken to be 10~ 15 and for a typical frequency range centered 
around the i^bbo — Hz. Furthermore, concerning the results of Figs. [2] and [3] the following 
comments are in order: 

• the pulsar bound, because of the smallness of its frequency range, is not directly 
relevant to constrain the secondary contribution: in Planck units the result would be 
A a < 10 12 for n a > 2 (see Fig. El plot at the left); 
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a ft 

Figure 2: The nucleosynthesis limits on the amplitude A a (expressed in units of the Planck 
mass Mp) for different values of n a and of /3. 

Q M0 =0.1 299, A R = 2.28x1 0" 9 , r T = 0.36, n s = 0.982 hjj Cl m =0.1299, A R = 2.28x 10~ 9 , r T = 0.36, n s = 0.982 
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Figure 3: The pulsar and Ligo/Virgo bounds as well as the expected sensitivities of planned 
space-borne interferometers such as Lisa and Bbo/Decigo. As in Fig. MA a is expressed in 
units of the Planck mass Mp. 

• the Ligo/Virgo bound is more relevant insofar as it implies A a < 10 5 ' 6 in Planck units 
for n a > 3: still since the frequency is too small and the corresponding window too nar- 
row the relevance of the obtained constraint is not comparable with the nucleosynthesis 
bound (see Figs. [2] and El plots at the left in both figures); 

• the Ligo/Virgo bound is moderately effective in constraining the maximal frequency 
of the secondary contribution for fixed values of n a : from Fig. [3] (plot at the right) 
(3 < 0(6) implies A a < C(10 5 ) for the case of violet spectral indices; 

• the big-bang nucleosynthesis bound, as expected, is the the most effective in constrain- 
ing both A a as well as (3 in the interesting ranges of spectral indices; in particular the 
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upper limit on A a is C(10 - 3 ) for n a > 2 (see Fig. [2]). 

By looking at Fig. [3] (plot at the left) it is possible to spot a corner of the parameter 
space where A a is compatible with the nucleosynthesis bound and where, simultaneously, 
the putative signal could be detected by the forthcoming space-borne interferometers. This 
corner of the parameter space is, roughly speaking, for 1 < n a < 2 and 1CT 2 < A a < 1. 
It should be however borne in mind that the results illustrated in Fig. [3] depend on a 
specific value of r T . Figures [2] and [3] present a scanning for different values of ft > 1. In the 
complementary case (i.e. ft < 1) the maximal frequency z/ max (see Eq. ( II. 3p ) gets smaller. 
In spite of the latter occurrence the results discussed so far are only slightly modified in the 
interval 10~ 5 < ft < 10. Indeed the big-bang nucleosynthesis bound is sensitive to all the 
frequencies of the spectrum for v > but a change in ft only changes z/ max which is the 
ultraviolet cut-off buts also the normalization frequency of the waterfall spectrum. As far 
as the Ligo/Virgo bound is concerned, its relevance becomes slightly more pronounced for 
ft < 1. For instance for 10~ 4 < ft < 10 -3 the Ligo/Virgo bound would demand A a < 10 3 
in Planck units and in the case n a = 4; for the same spectral index, Fig. [3] would imply 
A a < 10 6 always in Planck units. 

Before closing the present section few final remarks are in order. The second-order contri- 
bution of the metric fluctuations can also produce a secondary background of relic gravitons 
whose spectral slope and amplitudes are, however, fixed by the slope and amplitude of the 
curvature perturbations (i.e. n s and An) as discussed in [21]. Since n s < 1 the resulting 
secondary spectra will contribute to smaller frequencies and will represent a (subleading) 
correction to the quasi-flat plateau of the primary contribution to the spectral energy den- 
sity whose features will be briefly touched upon in section Therefore the second-order 
contribution of the metric fluctuations has been neglected both because it is proportional to 
*4jj and because it modifies preferentially the small frequencies. 

Equations (13. ip and (13. 2p do not account for the suppression of the quasi-flat plateau, at 
intermediate frequencies, due to the neutrino free streaming [32J. The evolution of the rela- 
tivistic degrees of freedom as well as the late dominance of dark energy have been neglected 
|33j . As it will be explicitly shown in the next section in the combined effect of neutrino 
free streaming, of the late dominance of dark energy and of the evolution of the number of 
relativistic species affects the amplitude of the intermediate frequency range by, roughly, one 
order of magnitude and this is the reason why these effects, even if conceptually important, 
are often ignored in explicit analyses (see, e.g. [H]), given our ignorance of the value of 7*t- 

4 High-frequency gravitons in the concordance model 

The primary and the secondary contributions to the spectral energy density of the relic 
gravitons will now be illustrated using the numerical techniques described in (3] (see also 
|10j ) and enforcing the constraints derived in section [3j To avoid confusions, when talking 
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about the concordance model we refer to the minimal extension of the ACDM paradigm which 
includes tensors. The latter extension is realized, as swiftly mentioned after Eq. (13.31) . when 
the spectral index does not run and rx = 16e fixes, at once, the amplitude and the spectral 
index of the tensor power spectrum. In Fig. [4] the primary and secondary contributions are 



£1 M0 = 0.1299, A R = 2.28 x 1CT 9 , r T = 0.36, n s =0.982 hjj £i M0 =0.1299, A R =2.28x 10~ 9 , r T = 0.1, n s = 0.982 




Figure 4: The primary and the secondary contributions to the spectral energy density are 
illustrated in the case of the parameters reported in Eq. ( I3.7P (plot at the left) and in the 
case of the parameters of Eq. (I3.7I) but taking into account also the neutrino free streaming 
and the late-time contribution of the dark energy. 

illustrated as they have been derived in the section [3J The maximal amplitude of A a has 
been taken to be 1 in Planck units. In the left plot of Fig. H] the effect of neutrino free 
streaming and the late-time dominance of the dark energy have been neglected. Conversely 
in the right plot both contributions have been included. 

The overall effect of the neutrino free streaming is to suppress the intermediate quasi-flat 
plateau for z/ eq < v < z/ bbn . Indeed, assuming that the only collisionless species in the thermal 
history of the Universe are the neutrinos, the amount of suppression can be parametrized by 
the function [32] (see also [15]) 

S(R U ) = 1 - 0.539^ + 0.134R* (4.1) 

where R v is the fraction of neutrinos in the radiation plasma, i.e. 

r / N \ 

Ru = —^, r = 0.681 (-^J, R 7 + R v = l. (4.2) 

In the case R v = (i.e. in the absence of collisionless patrticles) there is no suppression. 
If, on the contrary, R v ^ the suppression can even reach one order of magnitude. In the 
case N v — 3, R v — 0.405 and the suppression of the spectral energy density is proportional 
to S 2 (0.405) = 0.645. The suppression due to neutrino free streaming will be effective 
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for relatively small frequencies which are larger than v eq and smaller than the frequency 
corresponding to the Hubble radius at the time of big-bang nucleosynthesis, i.e. z/bbn- 

By carefully comparing the left and the right plots it is also possible to infer a global sup- 
pression which is actually due to the late-time dominance of dark energy. In the concordance 
model (i.e. the ACDM model) the dark energy is parametrized in terms of a cosmological 
constant and, consequently, the redshift of A-dominance will be 1 + z\ ~ (^/^mo) 1 ^ 3 - The 
mode reentering the Hubble radius at ta will have a frequency 

u A = 2.629 x ItW^U (%V /3 (4U 1/3 Hz. (4.3) 
V0.735/ V0.243/ V0.757/ v ; 

The adiabatic damping of the mode function across t\ reduces the amplitude of the spectral 

energy density by a factor (Q M0 /fl\) 2 ■ For the typical choice of parameters of Eq. (13. 7p 

the overall suppression of the spectral energy density is (9(0.103). There is, in principle, a 

further modification of the spectral slope occurring between 0.3 and 0.2 aHz (laHz = 10~ 18 

Hz) which has been taken into account but which is quantitatively immaterial. 

Finally the overall suppression of the curves of Fig. H] receives also a (subleading) contri- 
bution from the evolution of the number of (effective) relativistic degrees of freedom of the 
plasma entering the spectral energy density both as g p (i.e. the spin degrees of freedom per- 
taining to the total energy density) and as g s (i.e. the spin degrees of freedom pertaining to 
the total entropy density). In principle if a given mode k reenters the Hubble radius at a tem- 
perature Tfc the spectral energy density of the relic gravitons is (kinematically) suppressed 
by a factor which can be written as (see, for instance, [33]) (g p (Tk) / g p o)(g s (Tk) / g s o)~^ 3 ■ At 
the present time g p o = 3.36 and g s o = 3.90. The maximal expected suppression will be of 
the order of 0.38 in the standard model where g p = g s = 106.75 for temperatures larger than 
the top quark mass. In popular supersymmetric extensions of the standard model g p and g s 
can be as high as, approximately, 230. This will bring down the figure given above to 0.29. 

The spectra computed in the concordance model the last remark is that the quasi-flat 
plateau has been computed long ago by different authors [3U ESJ EE] and under different 
approximations. Thorough the years various completions of the standard paradigm have 
been proposed but none of them provides a strong gravitational wave background either at 
small or at large frequencies. An exception is represented by the presence of and early stiff 
phase [2"0"l I3"T] which can arise also in braneworld scenarios [3B] as well as in other contexts 
[39] . In Fig. [5] the primary and secondary contributions to the relic graviton spectra in the 
concordance model is compared with the spectral energy density computed in the case of the 
what has been dubbed TACDM model in [10]. The TACDM model (for tensor-ACDM) is an 
early time completion of the standard ACDM model where, after inflation, the dominance of 
radiation is delayed by a phase where the total speed of sound of the plasma is larger than 
the speed of sound of a radiation fluid. In terms of the barotropic index we shall have that, 
during the stiff phase, w t > 1/3. In Fig. [5] (plot at the left) the case w t — 1 is compared 
with the case typical of a waterfall field (i.e. n a = 4) and for the same range of parameters 
of Eq. ( 13 . 7j) . The case w t = 1 corresponds to a total sound speed coinciding with the speed 
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Figure 5: The secondary contribution to the spectral energy density of the relic gravitons 
is compared with the primary contribution computed in the TACDM scenario, i.e. when 
the standard radiation-dominated phase is preceded by a stiff epoch characterized by a total 
speed of sound larger than the one of a radiation-dominated plasma. 

of light and the spectral energy density of relic gravitons induced in such a case have been 
discussed in various reprises [20]. For 1/3 < w t < 1 the situation is partially illustrated the 
right plot of Fig. [5] where different values of rx are also allowed. 

5 Concluding remarks 

The secondary contribution to the spectral energy density of the relic gravitons induced by 
a waterfall-like field has been computed and constrained. The secondary contribution has 
been also combined with the standard signal arising in the case of the concordance ACDM 
paradigm and in its neighboring extensions such as the TACDM model. 
The obtained results suggest the following series of considerations: 

• for small frequencies (in the aHz range) the primary contribution of the concordance 
model remains unchanged by the secondary contribution; 

• for large frequencies, in the GHz range, the secondary contribution can be even 7 or 8 
orders of magnitude larger than the nominal amplitude obtained, in the concordance 
model, when the tensor spectral index is not allowed to run and the WMAP 7yr 
constraints on rx are implemented; 

• for intermediate frequencies, compatible with the observational window of wide-band 
interferometers the potential signal dies off faster than in the case of the TACDM 
paradigm where a broad high-frequency spike arises because, right after inflation, the 
speed of sound is stiffer than the speed of sound of the radiation plasma; 
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• the explicit constraints on the amplitude and on the spectral index of the waterfall 
field indicate that the nucleosynthesis constraint is more stringent than the direct 
Ligo/ Virgo bounds as well as more stringent than all the other bounds usually appli- 
cable to relic graviton backgrounds; 

• the secondary contribution is still ineffective for frequencies coinciding either with 
the (putative) Lisa and Bbo/Decigo sensitivities, however, if the spectral index is 
sufficiently small (i.e. the high-frequency hump sufficiently broad) also the secondary 
contribution can be relevant also at intermediate frequencies by so overlapping with 
some of the predictions of the TACDM model; this tuning is however not theoretically 
justified at present. 

The present results can be improved in a number of ways which have been partially mentioned 
in the bulk of the paper. In general terms steep spectra of waterfall-like fields might turn out 
to be a potential candidate for enhancing the spectral energy density of the relic gravitons at 
high frequencies. Along a more conservative perspective the obtained results will contribute 
to a better quantitative understanding of the parameters characterizing waterfall-like fields 
in conventional and unconventional hybrid models. It is interesting to remark that if the 
concordance model is complemented by a large tensor contribution at high-frequency low- 
frequency determinations of 7~t and high-frequency limits on the spectral energy density 
might not independent. In the future these kinds of models must be either ruled out (or 
ruled in) by combining low-frequency and high-frequency determinations of the spectral 
energy density. 
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